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Controllability of spin-boson systems 

Ugo Boscain, Paolo Mason, Gianluca Panati, and Mario Sigalotti 


Abstract 

In this paper we study the so-called spin-boson system, namely a two-level system 
in interaction with a distinguished mode of a quantized bosonic field. We give a brief 
description of the controlled Rabi and Jaynes-Cummings models and we discuss their 
appearance in the mathematics and physics literature. We then study the controllabil¬ 
ity of the Rabi model when the control is an external field acting on the bosonic part. 
Applying geometric control techniques to the Galerkin approximation and using per¬ 
turbation theory to guarantee non-resonance of the spectrum of the drift operator, we 
prove approximate controllability of the system, for almost every value of the interaction 
parameter. 


1 Introduction and description of the models 


1.1 The Rabi model alias the standard spin-boson model 

In the Hilbert space "H = h^(R, C) ® C^, we consider the Schrodinger equation 

i^tV^ = ilRabiV', (1) 


where 

h^Rabi = + X^)®t + -l®az + gX®ai, (2) 

for some hxed positive constants uj and P, and 5^ G R, where X is the multiplication operator 
dehned by {X'ip){x) = x'lplx), and P = —idx- Here and in the following we use cri,cT2,cr3 to 
denote the Pauli matrices 


CTi = 


0 1 
1 0 


0-2 = 


0 -i 
i 0 


= 


1 0 
0 -1 


Using the natural isomorphism Ti = h^(R, C^), one rewrites Q as the system 




+ x )^/J+ + -V'+ + 




idt'ip- = + x)i)-- -iIj- + gxi)+. 


In the physics literature, it is customary to introduce the raising and lowering matrices 


= 2(^1 + 1^2) 



= 2(^1 “ 1^2) 


0 0 \ 

1 oj- 


1 


Setting 


a — + i-P), 


at = ^(X-iP), 


and omitting tensors, Pnabi appears in physics literatnre as 

^^Rabi = ^ a^)(cr- + o-+)- 


( 3 ) 


Modeling of spin-boson interactions was initiated by Rabi in the 30 ’s im |22]. To our 
knowledge, the Hamiltonian (|^ was first derived from a more fundamental model in the 
milestone paper jl3j . where also the simpler Jaynes-Cummings Hamiltonian appears (see 


Section 1.2). A throughout discussion on various aspects of this model can be found in |3]. 
Recent results on its integrability have been obtained in j^. 

The Hamiltonian ([^ is ubiquitous in the literature, since it encodes one of the simplest 
possible couplings between a harmonic oscillator and a two-level system. The physical inter¬ 
pretation of the two factors appearing in the tensor product varies according to the context. 
We briefly outline the interpretation in cavity QED in the following digression. 


1.1.1 Digression: the Rabi model in the context of cavity QED 

While we refer to the specihc literature for an exhaustive treatement of the problem PI CUES], 
we briefly sketch the derivation of the Hamiltonian ^ in the context of Cavity Quantum 
Electro Dynamics (CQED). In a typical CQED experiment, an atom moves across a Fabry- 
Perot cavity interacting with the quantized electromagnetic (EM) held of the cavity. The 
atom is initially prepared in a special state, to guarantee both a strong interaction with the 
cavity held and a decay time (to the atomic ground state) longer than the experiment time- 
scale. The hrst goal is obtained by choosing a hydrogenoid atom (e.g., a Rubidium atom) 
prepared in a circular Rydberg state with high orbital angular momentum, and so with a large 
dipole moment. The second goal requires instead that transitions to the other atomic states 
are not resonant with the radiation from the environment, typically in the visible or infrared 
spectrum. Since the atomic energies are (where is fhe Rydberg constant, 

corrected for the Rubidium), one select n G N in such a way that the transition frequency 
f^at := — En-i) is iu the range of the microwaves (in the case of Rb, this is obtained by 

choosing n = 51). The cavity is designed in such a way that a normal mode of the EM held 
has frequency a;, which is almost resonant with flat, i-e-, |f^at — cal a;,f2at- 

Under these conditions, the only appreciable interactions are those among the two eigen¬ 
states |n) and |n — 1), and the distinguished mode of the EM held with frequency ca. The 
linear space generated by \n) and |n — 1) corresponds to the factor in EL, while the quan¬ 
tized EM mode yields a quantum harmonic oscillator, corresponding to the factor T^(R, C). 
This explains the structure of the hrst two terms in ([^. 

A deeper understanding of (|^ is obtained by deriving it, heuristically, from the Pauli-Fierz 
model, see |29l Chapter 13]. We consider a hydrogenoid atom, as, e.g.. Rubidium, interacting 
with the quantized electromagnetic (EM) held in the cavity. The electron, with ehective mass 
m* and charge —e, is described by a form factor (p G C'“(R^), which is assumed to be a radial 
positive function, normalized so that ip{x)dx = 1. The atomic core is supposed at rest at 
the origin or, equivalently, moving across the cavity with constant velocity. 
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Within the Pauli-Fierz model, the quantized EM fields are described by operator-valued 
distributions acting on the Hilbert space "Hat ® "Hf, where "Hat = L‘^(R,^,dr) corresponds to 
the atom, while 

nSN 

is the bosonic Fock space corresponding to the field j29]. In the Coulomb gauge, the vector 
potential (in free space) reads 

Mr) = E / M) exik) ® a(A;, A) + 0 a{k, A)t) , (4) 

^ Ar 3 y 


where uj{k) = c\k\, (p denotes the Fourier transform of 99 , the three vectors {k/\k\, ei(/c), e 2 {k)} 
form an orthonormal basis at every point k G R^, and a(fc. A) and a(/c, A)1 satisfy the canonical 
commutations relations. Analogously, the transverse electric field reads 

A'y,±(r) = E [ dk \l^^^-(p{k) ex{k) i (e'^’’’ ® a(/c. A) - ® a(fc, A)^) . (5) 


Actually, since the field is confined in the cavity, described by a compact set C C R^, the 
integral appearing in (|^ and ([^ can be replaced by an infinite sum where X labels the 

solutions to the eigenvalue equation AA = A^A satisfying V-A = 0 and appropriate boundary 
conditions on dC, see jSSl Sec. 13.6]. 

The (spinless) Pauli-Fierz Hamiltonian is the self-adjoint operator 

H = {ihS/r + + 1 ® Rf - hy,,Coui(r) ® 1, ( 6 ) 


where 


= ( dk hu{k) a{k, XYa{k, X) 

A=1 


corresponds to the field energy, while 


h^jCoul (^) 







ip{x-r)ip{y) 
47r|a; — y\ 


is the Coulomb interaction smeared through ip. 


First, since the electron is bound to the atomic core, “one loses little by evaluating the 
vector potential at the origin instead of at r, the position of the electron ” |29l Sec. 13.7]. 
Within such approximation, and after a suitable change of gauge, the Hamiltonian becomes 

Rdip = + Kflouiir)^ + Hi - er ■ E^,±(0) edip(r), (7) 

Obviously, in the comparison, the order of factors in the tensor product should be reversed, so that 
L2 (R,C)oC 2 is an effective reduced space for Hf OTfat- Since both conventions are well-established, we did 
not dare to reverse the order of factors in any of the two models, namely in ([^ and (|^. 


To derive the Rabi Hamiltonian (1^ from (1^, several approximations are in order 


( 1 ) 
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where r ■ -E(O) = ® The expression for hTdip involves the dipole operator 

D = —er, thus explaining the name dipole approximation. Here a quadratic potential appears, 
namely 

edip(r) = ^ j^dke^\ip{k)\^^ r^, 

which is usually neglected in this context. 

As a second crucial approximation, only a distinguished mode of the cavity field is con¬ 
sidered as interacting with the atom. If fc* is its wavevector, the integrals (or sums) above 
reduce to the evaluation at /c*. Accordingly, and assuming that (p{kf) = 1, one has 


^ y 2a;(fc ) [1 ® (a(fc*,A)-fa(/c*,A)'^)] 

eA(fcJ [l ®i(a(fc*,A)-a(A;*,A)^)] . (8) 

A=1 ' 

2 

yy fujj{k^)a{K,X)^a{K,X). 

A=1 

Neglecting the sum over the polarization vectors ca’s, we select an index A* G {1,2} and 
we set uj{k^) = u and a(/c*,A*) = Oc, where the subscript c stands for cavity. Then one has 
the correspondences 

A^(0) ~ Oc aj = V2Xc, and .^<^,±(0) ~ i(ac - of.) = -V2Pc, (9) 

where the last equality refers to the fact that the operators -^{cic + aj) and ^(cic — of) can 
be identihed, up to unitary equivalence, with a canonical Schrodinger pair (Xc, Pc) acting 
on L‘^(R,dq). Next, an inverse Fourier transform : L‘^(R,dq) —)■ L‘^(R,dx) is per¬ 

formed, which intertwines the pair {Xc,—Pc) with the pair (P, X), where P = —id^ and 
{X'il)){x) = x'ip{x). Correspondingly, OcP = —in, so that Pf is unitarily transformed into 
Pf = hwctt^a, which corresponds to the hrst addendum in (|^, up to an irrelevant constant. 
Moreover, 

A^piff) ~ —i(a — a^) = V^P, and P(^^j_(0) ~ a -|- = \/2X, (10) 

which shows that the interaction term (the third addendum in ([^ or in ([^) is proportional 
to the electric held, which is assumed to be uniform in view of the dipole approximation. 

Third, one assumes that effectively only two atomic states |n — 1) = \g) and \n) = |e) are 
appreciably coupled to the radiation held. Then, the factor Pat is replaced by a 2-dimensional 
Hilbert space 

PsL = Spauc {|^), |e)} ~ 

and the atomic Hamiltonian (the hrst term in ([^) by the 2-level Hamiltonian where 

hflat En Eji—i. 

Finally, if the two selected atomic states are circular Rydberg states one has that the dipole 
operator is oh-diagonal in the basis (l^f), |e)}, namely 

(^|T>|^) = 0, (e|P|e) = 0. 


^^( 0 ) = 
E^,±{0) = 

Pf 
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Moreover, there exist d > 0 and a versor Cat G such that de^t = (e| D \g). We identify a 
state V’at ^ H2L with the pair (-^g, ■^e) ^ via V’at = i^g \g) + V’e |e), so that 

= d (V’g |e) + Cat [s')) • 

By using the latter fact, one obtains 


D ■ ^y,,±(0) = ex,{K) ■ [e*t (a+0 i(ac - aj)) + Cat (o-_ ® i(ac - a];))] , (11) 

where, with a little abuse of notation, we wrote a- = \g) (e| and <7+ = |e) {g\. When adapting 
the previous formula to concrete physical situations, by specializing to the actual values of 
Cat and e\{K), one obtains some terms analogous to the ones appearing in ([^, possibly with 
different weights. In particular, if e*^ = Cat, one obtains a coupling term proportional to 
(Ji (S) Pc, which, after the Fourier transform mentioned above, becomes ui ® X and so agrees 
with the third addendum in (j^ or in ([^. 

This concludes the heuristic “derivation” of the Rabi Hamiltonian, starting from a fun¬ 
damental model as the one given by ([^. On the other hand, a rigorous derivation of the 
Hamiltonian (|^ from a mathematical model of non-relativistic QED, as, e.g., the Pauli-Fierz 
model, in a suitable scaling limit, is to our knowledge absent from the literature. 


1.2 The Jaynes—Cummings model 


The Rabi model presented in the previous section has been considered difficult to treat in 
view of the widespread opinion that it is not explicitly solvable (cf. [7]). 

For this reason Jaynes and Cummings jT3] proposed an approximated Hamiltonian, which 
is obtained from i^Rabi by neglecting the so-called counter-rotating terms. More precisely, the 
interaction term in iJRabi, namely 

(a^ -|- a)cri = (aV_ -|- aa+ + aV+ -|- acr_), 

is replaced by a new interaction term where the last two terms in the right-hand side are 
neglected under the assumption that 


|H - cal < fl. 

According to this procedure, one obtains 


g <^VL,oj. 


( 12 ) 


(13) 


//jo = w (ota + A + ^<^3 + A (“*''- + “'^+) 

= ^(/’^ + 2f^) 0 1 + 1 ® ^03 + |(A' ® <Ti - F 0 cTj). 

This model is considerably simpler than the Rabi one, since it admits the conserved quantity 

C' = a^a®l + 10( J Q 
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which represents the total number of excitations. Indeed, the interaction term in ifjc contains 
two parts: the term aV_ creates an excitation of the bosonic mode and destroys one of the 
two-level system, while acT_|_ acts in the opposite way. 


If the almost-resonance condition (12) is satisfied, the latter oscillations average to 

-1 


The heuristic justihcation of the approximation leading to (13) is based on separation of 
time scales. Indeed, when considering the dynamics in the interaction picture, the terms aV_ 
and 0(7+, which conserve the total number of excitations of the system, evolve periodically 
with frequency |f2 — a;|, while the remaining terms aV_|_ and acr_ oscillate with frequency 
If -|- OJ. 

zero on the longer time scale |f2 — a;| 

On the other hand, a rigorous mathematical justihcation for this approximation seems absent 
from the literature. However, as already noticed by Rouchon in |27], it seems possible to 
obtain it by adapting the methods developed in pUl IT^ I3U] . A similar task, by analytical 
methods, has been accomplished by Ervedoza and Fuel on a related model |12j . 

Notice that in applications to circuits QED the assumption 0,0; is in general not 
satished, since it is possible to achieve spin-hied interactions which are not much smaller than 
the held energy j28]. This is a motivation for the direct study of the Rabi model, which is 
pursued here. 


1.3 Modeling an external control field 

In a wide variety of experimental situations one can act on the system by an external held. 
The goal of the controller could be to lead the system from a given initial state to a prescribed 
hnal one. For spin-boson models, this amounts to study the control problem 


(14) 

where Hq represents either the Rabi or the Jaynes-Cummings Hamiltonian and H^. is a self- 
adjoint operator describing the coupling between the system and the controlled external held. 
The operator depends on the control function u which, in general, takes values in C™. 

In most cases the external held can act on the bosonic mode only, while the spin mode is 
not directly accessible. This leads to a control Hamiltonian of the form 


Hc{u{t)) = hc{u{t)) 0 1. 


(15) 


where hc{u{t)) is a self-adjoint operator in L^(R). One of the simplest form for the operator 
hc{u{t)) is the following 

hc{u{t)) = u{t)X. (16) 


As argued in Section 1.1.1, in the context of cavity QED and within the dipole approxi¬ 


mation, the operator X corresponds to the electric held at the origin, see (10). Therefore, in 


this context, the control (16) corresponds to an external EM held which rescales the value of 


E^ _i_(0) by a time-dependent factor u{t). Notice that the linearity in X is a consequence of 
the dipole approximation. 

Simmetrically, one might also consider a control consisting of an external held which mod¬ 


ulates the value of A^{0), as e.g., an externally-controlled magnetic held. In view of (10), this 
amounts to consider a control term hc{u{t)) = u{t)P, a problem which will be investigated 
elsewhere. 
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1.4 Other related controlled models 

Controllability of finite-dimensional approximations of spin-boson systems have first been 
obtained in the physics literature, via constructive methods, in d! and jT^. The results have 
been extended to the full infinite-dimensional model in [3l] via analysis of non-resonances and 
in jl] by adiabatic arguments (see also i). 

A mathematical description of the controllability properties of a related model is given in 
|12) . The main technical tool is an explicit estimate of the approximation error with respect 
to the Lamb-Dicke limit. 

For the counterpart of the Eberly-Law model for more than one trapped ion an approxi¬ 
mate controllability result is obtained in |1] (see also 1211 ), based on the analysis on a sequence 
of nested finite-dimensional systems. In the recent papers [T^ [T8] the (approximate) control¬ 
lability of the system is established by considering different families of controlled dynamics. 


2 Approximate controllability of the controlled Rabi model 


We consider here the approximate controllability problem for a system of the form (14), where 
Ho is the Rabi Hamiltonian and takes the form ([l^-([l^. 

The problem under consideration is then 


{ idtiJ = + u{X 0 1)V’, 

u e [0,5], 

-cul -C fl, 

5, > 0. 


(17) 


The main result of the paper is the following. The precise definition of approximate 
controllability will be given in the next section. It basically means that, for every choice of 
the initial and final state, there exists an admissible control law u, depending on the time, 
which steers the initial state arbitrarily close to the hnal one. 


Theorem 1. System (17) is approximately controllable for almost every g eH. 


The proof of Theorem goes by studying the applicability of a general approximate con¬ 
trollability result in dependence on the parameter g. In order to do so, we have to use 
perturbation theory in the parameter g up to order 4. 


3 An approximate controllability result 

We are going to recall a general controllability result for bilinear quantum systems in an 
abstract setting. 

In a separable Hilbert space H, endowed with the Hermitian product (•,•), we consider 
the following control system 


d 

dt 


= {A + u(t)B)'ip, 


u(t) e U, 


where {A, B, U) satisfies the following assumption. 


(18) 
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d)—( p &—0- 

Figure 1: Each vertex of the graph represents an eigenstate of A (when the spectrum is not simple, 
several nodes may be attached to the same eigenvalue). An edge links two vertices if and only if B 
connects the corresponding eigenstates. In this example, {4>i,B(j>2) and {(f>i,B(p 3 ) are not zero, while 
{(j)l,B(j)4) = Bcps) = {4>2, B(j)4) = 0. 


Assumption 1 . U is a subset of R and {A,B) is a pair of (possibly unbounded) linear 
operators in H such that 

1. A is skew-adjoint on its domain D{A); 

2. there exists a Hilbert basis (0fc)fceN of H consisting of eigenvectors of A: for every k, 
A(j)k = with Xk in R; 

3. for every j in N, cfj is in the domain D{B) of B; 

4- A + uB is essentially skew-adjoint for every u G U ; 


5. {B(j)j,(j)k) = 0 for every j, k in N such that Xj = Xk and j ^ k. 

If {A,B, U) satisfies Assumption then A + uB generates a unitary group t i—)■ _ 

By concatenation, one can define the solution of (18) for every piecewise constant function u 
taking values in U, for every initial condition ^jq given at time to. We denote this solution by 


For every j,k G N, we denote bjk = {4>j,B(j)k)- A pair {j,k) in is a non-resonant 
transition of (A, B) if bjk ^ 0 and, for every I, m, |Aj — Afc| = |Ai — Am| implies {j, k} = {I, m} 
or {/, m} n {j, k} = 0. 

A subset S of is a chain of connectedness of (A, B) if for every j, k in N, there exists a 
finite sequence pi = j,p 2 ,... ,Pr = k for which {puPi+i) G S for every I and 7 ^ 0 

for every / = 1,..., r — 1. A chain of connectedness S of (A, B) is non-resonant if every (j, k) 
in S' is a non-resonant transition of (A, B). 


Definition 2. Let {A,B,U) satisfy Assumption We say that (18) is approximately con¬ 
trollable if for every e > 0, for every ^ H, there exists a piecewise constant function 

Ue : [OjTg] —)■ U such that ||T^ — f/’ill < 


Theorem 3 ([5]). Assume that [0,5] C U for some 5 > 0 and let {A,B,U) satisfy Assump¬ 
tion and admit a non-resonant chain of connectedness. Then system (18) is approximately 
controllable. 









4 Proof of Theorem [U 


In this section, we consider the approximate controllability problem for (17) and we prove 


Theorem The proof of the theorem is based on a suitable application of Theorem 

The strategy of the proof is the following. We first show in Subsection |4.1.1| that, for 
almost every g in R, some relevant pairs of eigenvalues of RRabi satisfy the non-resonance 
condition, see (22). This goal is achieved by exploiting the analyticity of the eigenvalues 
and by using perturbation theory. Then, in Subsection |4.1.2[ we prove that these pairs of 


eigenvalues correspond to non-resonant transitions, according to the definition above. The 


degenerate case a; = hi is treated separately in Section 4.2 


Preliminarily, we introduce some additional notations. Denote by i7Rabi,o the Hamiltonian 
-f^Rabi where we set g = 0. The general form of Ruabi is then 


-f^Rabi — -ffRabi,0 + QV, 


(19) 


where V is the operator X ® ai. 

Let be the standard Hilbert basis of L^(R, C) given by real eigenfunctions of 

-dl + x^, so that + = {2j + l)gjj and j^xg)j{x)g)j+i{x)dx = \/{j + l)/2 for j > 0. 

Let z/i = (1, 0)^, z/_i = (0,1)^ be the canonical orthonormal basis of C^. Based on ((pj)jgN, 
we obtain a Hilbert basis of factorized eigenstates j G N, s G { — 1,1}, of RRabi,o 

whose corresponding eigenvalues are 


-3,s - ( j + 2 ) + 2^‘ 


( 20 ) 


If D is not an integer multiple of oo then each eigenvalue Ej^s is simple. (See Figure]^) In the 
following, for ease of notations, we write in bold the elements of N x {—1,1}, and for every 
j G N X { — 1,1} we define n(j),s(j) in such a way that j = (n(j),s(j)). 

For g E H, denote by E?, j G N x { — 1,1}, the eigenvalues of Ruabi repeated according 
to their multiplicities, and by 4)|,j G N x { — 1,1}, an orthonormal basis of corresponding 
eigenstates. By a consequence of perturbation theory, detailed in the Appendix, we can assume 
that, up to the choice of a suitable labelling, each map H, 3 g t-E- {E?, $?) G R x (R,C) 


is analytic, with E^ = E^, where E^ is defined in (20). When E^ is simple we can assume, in 
addition, that = <I)j. 

Under the assumption ID —a;| <C D, the only case in which i7Rabi,o has multiple eigenvalues 
is when a; = D. If this is the case all eigenvalues different from the lowest one are double. 


Equation (32) allows to identify the splitting of the two-dimensional eigenspaces induced by 


the perturbation g i—)■ FfRabi.o + gV, as detailed in Section 4.2 


In order to study the first and higher-order derivatives of g ^ Ej^ at g = 0, it is useful to 
introduce the quantities 


hij (*hi, {x ® j ^n(i),n(j)-l 


^(j) 


+ h 




n(j) + 1 


(1 - ( 21 ) 
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4.1 Case cj 7 ^ Q 

4.1.1 Step I: Relevant eigenvalue pairs are non-resonant. 

Let us first prove that for almost every (7 G R and every i,j,k, 1 G N x { — 1,1}, with (i, j) 7 ^ 
(k, 1) and i 7 ^ j, one has E? — ^ E^ — Ef. In order to do so, we observe that it is enough 

to show that for fixed i,j,k, 1 G N x { — 1,1} as before, the set 

Sij,k.i = {9 I E! -EfytE^- B?} (22) 

is of full measure. By the analytic dependence on g of the eigenvalues of Rnabi, this is 
equivalent to say that g ^ Ef — E? and g ^ E^ — Ef have different Taylor expansions at 
^ = 0 . 

Let us consider the Taylor expansion 


E! = Ei + J2 g-El 


(m) 


m=l 


The computation of the coefficients carried on below is based on the Rayleigh-Schrodinger 
series (see, for instance, [26l Chapter XII]). 

First of all we observe that Ei — Ej = E]f_ — Ei is equivalent to n(i) — n(j) = n(k) — n(l) 
and s(i) — s(j) = s(k) — s(l), in the case in which 12 is not an integer multiple of a;/2. If 
12 = (2m + l)ci;/2 for some non-negative integer m, then Ei — Ej = Ei^ — Ei implies 

Orn - 1 — 1 

n{i) + n(l) - n(j) - n(k) = —-— (s(j) + s(k) - s(i) - s(l)), 


and thus, if the left-hand side is an integer number different from zero, it must be |s(j) -|- 
s(k) — s(i) — s(l)| = 4, that is s(j) = s(k) = —s(i) = —s(l). 


The term E^^'^ coincides with Vj j 
= 0 for every j. 

Following [2n] we have that 


(<hj, (X ® (Ji)<hj), thus we deduce from (21) that 




Thus 




( 2 ) 


5^(F;^-F;j)-i(l-5,(j),q^))Mh 

V 


n(j),n(m)-l 


n(j) + 1 


+ 


n(j),n(m)+l 


^(j) 


- Ei)-' - Bj)-' ^ 

-{u: - 

u + g(j)12(2n(j) -11) 

2 ( 122 -n; 2 ) 
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Notice that the computation above is correct also for n(j) = 0, even if in this case 
is not a well defined eigenvalue of -ffRabi,o- Indeed, in this case the term 


-^j) 




counts as zero. 

(2) (2) (2) (2) 

Let us identify the values i, j, k, 1 such that under the assumption 

that Ei — Ej = E]f_ — El. Recall that we also assume that (i, j) ^ (k, 1) and i 7 ^ j. From the 

( 2 ) 

above expression of E^ ^ we have 

s(i) (2u(i) + 1) - s(j) (2n(j) + 1) = s(k) (2u(k) + 1) - s(l) (2n(l) + 1). (23) 

If n = (2m + l)a;/2 for some non-negative integer m and s(j) = s(k) = —s(i) = —s(l) then 


(23) gives u(i)-|-u(j)-|-n(k)-|-n(l)-|-2 = 0, which is impossible being the addends non-negative. 

The remaining case is when n(i) — u(j) = n(k) — n(l) and s(i) — s(j) = s(k) — s(l), in 
which case 


s(i)n(i) - s(j)u(j) = s(k)n(k) - s(l)n(l). 


(24) 


Then, either s(i) = s(j), which implies s(k) = s(l) and then, by (24), s(i) = s(j) = s(k) = s(l). 


or s(i) = —s(j), which implies s(k) = —s(l) and then, by (24), n(i) -|- n(j) = n(k) -|- n(l). In 


the latter case it must be i = k and j = 1, which is excluded by assumption. Therefore the 
nontrivial quadruples i,j,k,l satisfying both the equalities Ei — Ej = — Ei and = 


42 ) 


(‘ 2 '] 

— El are those for which s(i) = s(j) = s(k) = s(l) and u(i) — u(j) = n(k) — n(l). 


Let us now evaluate the terms Fj ^ as in 


We have 


d = E - Bj)-Vj,„v„,pVpj - 

Since 14,b 7 ^ 0 only if s(a) = —s(b) it turns out that Vj,„i and I4i,p are different from 0 only 

/o\ 

if s(j) = s(p) = —s(m), but then Vp,j = 0. Thus, recalling that Vj,j = 0, we have F^j = 0 for 
every j G N x { — 1,1}. 


We are going to complete the proof that the set defined as in (22) has full measure by 

showing that if i, j, k, 1 G N x { — 1 , 1 } are such that (i, j) 7 ^ (k, 1), i 7 ^ j, n(i)—n(j) = n(k)—n(l) 
(which follows from Ei — Ej = i?k — Ei) and s(i) = s(j) = s(k) = s(l) (which follows from 


E} 


( 2 ) 


E{2) ^ ^(2) 

J ^ 


Ep^), then E^^'^ - E^^'^ ^ E^^'^ - E^^l 


44) 


n(4) 


The general formula for E^ 


( 4 ) 


,see 


IS 




m4j,P4j.q4j 

E rijU„U.,pVpj|(£™ - £j)-‘(£p - + (£„. - - £j)-‘l 

^ - Ei)-^E^ - £'j)-‘U„r„jb,pV'pj - 
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Since Vjj = 0, only the first and third term of the right-hand side must be evaluated. 

Let us compute the first term in ( [25| . In order to avoid null terms we must assume 
s(j) = —s(m) = s(p) = —s(q) and thus n(j) 7 ^ n(p). Therefore the only nonzero terms in the 
sum are given by n(j) = n(m) -|- 1 = n(p) + 2 = n(q) -j- 1 (if n(j) > 1 ) and n(j) = n(m) — 1 = 
n(p) — 2 = ? 7 ,(q) — 1. We have 


2 / V 2 
-1 f ’^(j) + ^\ f ^(j) + 2 


(- 0 J - sij)Sir‘(-2ojr‘(-oj - s(j)SJ)-‘ ^ 

+ (a;-s(j)fi)-^(2a;)-^(n;-s(j)f2) , 2 /V 2 

Notice that the formula is correct also in the case where n(j) = 0 or n(j) = 1. 

Let us now compute the third term in (25). As before, to avoid null terms we assume 
s(j) = —s(m) = —s(p). The nonzero terms in the sum are given by n(m) = n(j) ± 1 and 
n(p) = ? 7 ,(j) ± 1 thus we have to sum four terms. We have 

E = (-U, - s(j)!J)-7-c - s(j)n)-‘ 


(-u - s(j)Q) ^(cn - s(j)f 2 ) 
-f (w - s(j)f 2 )"^(-(:n - s(j)f 2 ) 


-1 ^ ^(j) + 1 
2 2 
-1 ^(j) ^(j) + 1 


(ca - s(j)f 2 ) Aa;-s(j)f 2 ) 


-1 f ’^(j) + 1 


By summing up all the terms one sees that, for fixed s = s(j), the term depends 
quadratically on n(j), i.e. 

= C'o(s(j)) + C'i(s(j))n(j) + C' 2 (s(j))n(j)^ 
where the coefficient C 2 (s(j)) is given by 


C' 2 (s(j)) = s(j) 


.,fi(n ;2 + 3ff2) 


2 (a;2 _ (22)3 


^ 0 . 


So, if i, j, k, 1 are such that s(i) = s(j) = s(k) = s(l) = s and n(i) — n(j) = n(k) — n(l), we 
have 


_ ^(4) ^ ^(4) _ ^( 4 ) ^ C'i(s)(n(i) - n(j)) + C 2 (s)(n(i )2 - ^(j)^) = 

Ci(s)(n(k) - n(l)) + C 2 (s)(n(ky - nilf) 
^ C2(s)(n(i)^ - nOf) = C2(s)(n(k)^ - nilf) 
C 2 (s)(n(i) n(j)) = C 2 (s)(n(k) n(l)) 

<¥=¥ n(i) = n(k) and n(j) = n(l). 


This concludes the proof that for almost every g G R and every i,j,k, 1 G N x {- 1 , 1 }, 
with (i, j) 7 ^ (k, 1) and i 7 ^ j, one has E? — E? 7 ^ — Ef. 
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Figure 2: The dashed lines connect eigenvalues of -ffRabi that are “coupled” by the controlled 
Hamiltonian when g = 0, while the dotted lines connect eigenvalues that are coupled by the 
controlled Hamiltonian for almost all 5 ^ 7 ^ 0. 


4.1.2 Step 2: Coupling of the relevant energy levels. 

The proof of Theoremis then concluded, thanks to Theoremif we show that the controlled 
Hamiltonian a; ® 1 couples, directly or indirectly, all the energy levels for almost all g 

More precisely, we show below that ($|, {x ® 7^ 0 for almost every G R for all j, k 

such that s(j) = s(k) and |n(j) — n(k)| = 1 or s(j) = —s(k) and n(j) = n(k). See Figure 
As before, it is enough to show that the corresponding Taylor series in g is nonzero. 

Set $1 = <Fj + Em=i We have 


^n(j),n(k)-l Y 2 ^ Y 2 j ' 

This is enough to say that (<F|, {x ® 1)^*^) 7 ^ 0 for almost every g for all j,k such that 
■sQ) = s(k) and |n(j) — n(k)| = 1 . 

The term can be characterized through the relation 


(^Rabi,o + gx^ cri)($j + ) + o{g)) = {E; + gE\^> + o((7))($j + + o{g)). ( 26 ) 


( 1 ) 




Regrouping the first-order terms in (26) we get 

+ (x (g) ui) = 0. 




Rabi.O ■ 


(27) 


Denote by H the orthogonal projection on the orthogonal complement to <Fj. Applying H to 

(RRabi,o - + n(a; ® (Ti)<Fj = 0. 


(| 2 ^, we get 


Notice that the orthogonal complement to <Fj is an invariant space for the operator Ruabi.o ~ 
£jl, which is invertible when restricted to it. We write (RRabi.o ~ fo denote its inverse 
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(whose values are in the orthogonal complement to <hj). Thus, 


) = -(i/Rabi,o - Eil)-^U{x ® ai)<hj = (a: ® 

¥j 

The linear term in the Taylor expansion of {$?, {x ® 1)*^^) with respect to g is given by 
(<hj, (x (g) + (<hf \ (x (g) l)<hk) = - ^i)“^(<hi, (x (g) cri)<hk)(<hj, (x g) l)<hi)+ 

¥k 

^(Ej - Ei)-^(^i, (x 0 ai)<hj)(<hi, (x 0 l)$k). 

¥j 


Taking into account only the nonzero terms in the first sum gives s(j) = s(l) = —s(k), 
|?7,(j) — n(l)| = 1 and |u(k) — n(l)| = 1, so for fixed j, k we only have (at most) two terms. We 
assume n(j) = n(k), thus we have 




(-U + s(k)Sl)-‘dtW + (IJ + s(k)SJ)-'^ 

CO + s(k)r2(l + 2?7,(k)) 

2(f22_a;2) • 


Similarly, 


Hence, 


($f\(a;gl)<hk) 


CO — s(k)r2(l + 2?7,(k)) 
2(f22 _a;2) 


\ (x 0 l)<hk) + (<hj, (x 0 = ^2 ^^2 


which is different from zero for every k. Thus the Taylor series of {<hj, (x ® 1)*^*^) is nonzero, 
which concludes the proof of the theorem in the case u ^ Q. 


4.2 Case oj = Q 

Let us assume in this section that ca = hi. Hence 

Ej,i = -E'i+1,-1 = + 1) 


for every j G N and = 0. 

Recall that, as g varies in R, a complete set of eigenpairs for i^Rabi can be given in the 
form (T;?,<h?)jgNx{-i,i} with each E?^ and depending analytically on g and = Ey Notice 
that each <hj can be taken in T^(R, R) g and that we can take = $ 0,-1 = ^o0 ^-i- 
The eigenfunctions are characterized by the relation 




(This follows by standard perturbation relations recalled in the Appendix, see equation (32).) 
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Since, moreover, and d’j+i _i are orthogonal and of norm one, one gets, up to signs, 

-1 “ *^i+i-1 


Let us fix, by convention, 

+ *^i+i-1 


72 


= 


72 


*0 

^i+i-1 


72 


1 - $ 


i+i-1 


72 


With this choice (32) yields 


7_ /?9 1 = /i +1 

dg ^’^' 9=0 


—E^ 1 

dg 9+l--ll9=0 


J + 1 


Assume now that 


77'U rpv r^U r^U 

— I (B? -£)?! = —I iE?-E?) 

dg>s=<>'^ Ji n’ dg‘ii-«' ■“ 


(28) 

(29) 


and let us prove that ji = js and j 2 = j 4 . Write = uki for / = 1,...,4 for some 
fci, ..., ^4 G N. Then (28) implies that ki = k 2 + p and k^ = k^ + p for some p G Z. 

Assume that fci,. .. , /C 4 7 ^ 0 (the other cases being similar). Equality (29) then yields 

Si 7^2 + P - S2-s/k2 = S3\/k4+ p - 54-7^ (30) 


for some Si,... ,54 G { — 1,1} and we are left to prove that si = S3, S 2 = S4 and ^2 = ^4. 
Without loss of generality we can assume that k 2 + p = max{A; 2 , ^2 +^,^(4,^4 + p} and that 


can then rewrite (30) as 


Si = 1. Then p and the left-hand side in (30) are nonnegative, which implies that S 3 = 1. We 


7^2 +P - S 2 \/h = -s/ki+p - S 4 \/kl 


(31) 


If S 2 = S 4 then, by monotonicity of fc h-)■ ^k + p ± \/k we deduce that ^2 = k 4 and we 
are done. Finally, if S 2 = —S 4 then either S 2 = —1, and then k^ > k 2 , or S 2 = 1 and 
then 7^7+P = 7^4 + P + 7^ + 7^ > a/P + 7fe > V ^2 + P, leading in both cases to a 
contradiction. 

In order to apply Theorem we are left to prove that the controlled Hamiltonian X (g) 1 
couples, directly or indirectly, all the elements of the basis (<h?)jgNx{-i,i}. 

First notice that 


(K?,!. (X ® 1)>I>"+,,,) 



jeN, 


meaning that X 0 1 couples all basis elements of the type $ 7 - 

Similarly, one has that ® ^ ^ every j G N. Finally, one easily 

checks that ($7’ ® ^ ^ every j G N, completing the proof. 
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Appendix: global real analyticity of eigenpairs 

In this section we detail the proof of the analyticity of the eigenpairs {E?, <hj) as functions 
of 5 ^ G R. This result is somehow folklore and follows from well-known theorems in linear 
perturbation theory (na l26]). Since the latter are spread over the literature, we prefer to 
sketch here the main ideas of the argument, which is also useful in other contexts, as e.g. for 
the degenerate Gell-Man and Low theorem (ilElIlo]). 

Recall that, given two self-adjoint operators Hq, IT on a separable Hilbert space, W is said 
to be Kato-small with respect to Hq if D{Hq) C D{W) and for every a > 0 there exists b > 0 
such that ||IT'0|| < a||iLo'0ll + for every G D{Ho). According to this definition, it is 
easy to check that V, defined as in ( |I^ , is Kato-small with respect to ifRabi.o- 

Proposition 4. Let Hq^W he self-adjoint operators on a separable Hilbert space H. As¬ 
sume that W is Kato-small with respect to Ho and that Hq has compact resolvent. Then 
Hg = Ho + gW is self-adjoint with compact resolvent for every g & H, hence it admits a 
complete orthonormal system of eigenfunctions {'f’j{g)}j£-N ond the corresponding eigenvalues 
{£'j(5f)}jgN have finite multiplicity and do not accumulate at any finite point. Moreover, up 
to a suitable choice of labelling, for each j G N the function g i—)■ {Ej{g),‘ipj{g)) is analytic 
from R to'R.x'H. 

Proof. Self-adjointness of Hg follows from [251 Theorem X.12]. Let us first prove that Hg has 
compact resolvent for every S' G R (we provide a direct proof; see also [TTl Theorem VIL2.4]). 
SincJ^ 

{Hg - il)-i = {Ho - il)-' {1-g W{Hg - il)-i) 

and {Hq — il)“^ is compact, we are left to prove that W{Hg — il)“^ is bounded. Notice 
that for every a > 0 there exists 6 > 0 such that ||IT'i/'ll < ^\\{Hg — il)'0|| -l- &||'0|| for every 
G D{Ho) = D{Hg — il). Setting = {Hg — il)“^0 for 0 G H, we get 

\\W{Hg - il)-V|| < a||0|| + b\\{Hg - il)-V|| < (a + 6)||</.||, 

where the last inequality follows from \\{Hg — il)r 7 ||||r 7 || > \{{Hg — il)g,g)\ > ||?7|p, for g G 
D{Hg -il). 

If Ej{g) is simple then, by Rellich’s theorem [26l Theorem XIL8], up to a reordering 
of the spectrum of Hgi for g' in a neighborhood of g, the eigenpair parameterization g' t-)- 
{Ej{g'),fjj{g')) is analytic near g. The previous result can be generalized to the case of 
eigenvalues of multiplicity m > 2. Indeed, if Ej^{g) = ■ ■ ■ = Ej^{g) and ji,... ,jm are distinct, 
we can assume, up to relabelling, that g' h -)■ {Ej.^{g'),'ipj.^{g')),... ,g' {Ej^{g'),f)j^{g')) are 

analytic near g ( |261 Theorem XIL13]). Moreover, 

{'^jii9),Wi)j^{g)) = dik-^EjXg), l,k = l,...,m. (32) 

(See, for instance, [2].) 

In order to describe the spectrum by a countable family of functions which are globally 
analytic on R, we should ensure that all such locally analytic functions Ej{-) can be extended 
indehnitely. 

Notice that {Hg — il)“^ maps % into D{Hg) C D{W), so that W{Hg — il)“^ is well-defined on the whole 
Hilbert space H. 
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Fix j and let / be the maximal interval containing 0 such that, up to relabelling, Ej is 
analytic on I. Let us prove that, for every i G (0, +oo), [—£,£] C I. 

Set a = Take 6 > 0 such that ||hFt/’|| < a||iLot/’|| + guv'll for every G D{Hq). 

Hence, for -0 G D{Hq) and g G 

\\W^P\\ < a\\iHo + gW)^p\\+a\g\\\W^P\\+b\\^p\\ 

and therefore 


IlfT^II <2a||i7,^||+2611^11, ge [-£,£]. 


It follows from (32) that, for every G / fl [—£,£], 


d 

dg 


Ej{g) 


\{'ipj{g),W'il:j{g))\ < ||fF^/'j(^)|| < 2a\\Hg'ipj{g)\\ + 2b\\'ipj{g)\\. 


Hence, 




— 2a|-Fj(5')| + 2b, 


for almost every G / fl [—£, £]. 


A standard application of Gronwall inequality yields that Ej{-) is Lipschitz on / fl [—£,£]. 
In particular the limits of Ej{-) at the boundary of / fl [—£,£] are well dehned and they are 
eigenvalues of the corresponding perturbed operators (see e.g. |26l Theorem XH.7]). Therefore 
Ej{-) can be extended indehnitely on [—£,£] and thus also on the whole real line. □ 
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